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Abstract. We discuss the role of boundary conditions in determining the physical 
content of the solutions of the Schrodinger equation. We study the standing-wave, the 
"in," the "out," and the purely outgoing boundary conditions. As well, we rephrase 
Feynman's +ie prescription as a time-asymmetric, causal boundary condition, and dis- 
cuss the connection of Feynman's +ie prescription with the arrow of time of Quantum 
Electrodynamics. A parallel of this arrow of time with that of Classical Electrody- 
namics is made. We conclude that in general, the time evolution of a closed quantum 
system has indeed an arrow of time built into the propagators. 

1 Introduction 

In physics, dynamical equations often have a differential form and are solved 
under various boundary conditions. This is also the case in Quantum Mechanics, 
whose dynamics is encoded by the Schrodinger equation. The purpose of this 
contribution is to discuss, in a somewhat sketchy way, how boundary conditions 
imposed upon the Schrodinger equation determine the physical content of its 
solutions. And vice versa: in order to obtain the solutions of the Schrodinger 
equation that describe a given physical situation, boundary conditions that fit 
the physical situation must be imposed upon the Schrodinger equation. 

In the non-relativistic domain, the time-independent Schrodinger equation 
is realized, in the position representation, as a second-order differential equa- 
tion. If the potential is simple enough, we can exactly solve the differential 
equation by means of the Sturm-Liouville theory This theory yields an 

eigenfunction expansion, a unitary operator that diagonalizes the Hamiltonian, 
and a direct integral decomposition. In its turn, the direct integral decompo- 
sition yields, along with some physical requirements, a Rigged Hilbert Space 
(RHS) . The eigenfunction expansion, the unitary operator that diagonalizes the 
Hamiltonian, the direct integral decomposition, and the RHS contain much of 
the spectral and physical informations of the Hamiltonian. For convenience, we 
shall refer to them as the RHS properties (abbreviated RSP) associated to the 
Hamiltonian 0]. At first sight, it may seem that the Schrodinger equation cor- 
responding to a given Hamiltonian generates just one single RSP. However, this 
is not necessarily so. There are certain types of boundary conditions that, when 
imposed upon the Schrodinger equation, yield different RSPs. To be more pre- 
cise, the standing- wave boundary condition (to be defined below), and the "in" 
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and "out" boundary conditions of the Lippmann-Schwinger equation generate 
three different RSPs. 

Tin- Gamow vectors are the state vectors of resonances |5l(il7l<Sli)ll()ll| . Like 
the standing-wave and the Lippmann-Schwinger eigensolutions, they solve the 
Schrodinger equation. At infinity, however, the Gamow eigcnfunctions satisfy 
the purely outgoing boundary condition. This condition selects the (complex) 
resonance spectrum of the Schrodinger equation. 

The time-dependent Schrodinger equation is time symmetric, the reason for 
which it is generally believed that Quantum Mechanics is time symmetric 
And yet the time evolution of individual atoms or subatomic particles seems 
to have some directness. For example, imagine that we want to compute the 
probability for a particle to go from an initial space-time location (x, t) to a 
final space-time location (x', t'). Our basic notions of causality dictate that this 
probability be zero when t' < t. However, as far as the Schrodinger equation is 
concerned, this probability is also non-zero when t' < t. In order to obtain the 
causal result, we have to use the retarded propagator G + (x, t; x', t') (see for ex- 
ample This retarded propagator automatically yields a causal probability, 
because it vanishes when t' < t: 

G+(x,i;x',t') = 0, t'<t. (1) 

Therefore, the time evolution given by G + (x, t; x', t') has an arrow of time: a 
particle travels forward in time, never backward. As is well known, the time 
evolution given by the (e.g., retarded) propagator is equivalent to the time- 
dependent Schrodinger equation subject to an (e.g., retarded) causal boundary 
condition. Hence, the arrow of time built into G + (x, t; x', t 1 ) stems ultimately 
from causal boundary conditions. 

Quantum Electrodynamics (QED) provides a glaring example of how the 
propagators carry an arrow of time. In QED, the Feynman propagator is con- 
structed by imposing Feynman's +ie prescription upon the time evolution of 
particles and antiparticles |T3 . Particles travel forward in time, whereas an- 
tiparticles "travel backward" in time. Clearly, this prescription builds an arrow 
of time into the Feynman propagator. But note that this arrow of time is intro- 
duced by means of Feynman's +ie prescription, which is a boundary condition. 
Thus, the arrow of time of the Feynman propagator also stems ultimately from 
causal boundary conditions. 

The organization of this contribution is as follows: 

(z) In order to obtain the standing- wave eigensolution (r\E) , the "in" Lippmann- 
Schwinger eigensolution (r\E + ), and the "out" Lippmann-Schwinger eigensolu- 
tion (r\E~), we shall impose upon the Schrodinger equation the standing-wave, 
the "in," and the "out" boundary condition, respectively. Each of the eigensolu- 
tions (r\E), (r\E + ), (r\E~) yields an RSP of its own. Hence, we shall conclude 
that each of the eigensolutions (r\E), (r\E + ), (r\E~) has a physical content of 
its own. 
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(it) It will be apparent that what makes (r\E), (r\E + ), (r\E~) different from 
each other is the boundary conditions that they satisfy. More precisely, what 
makes them different from each other is their asymptotic behavior at infinity. 

(Hi) We will relate the asymptotic behavior at infinity of (r\E), (r\E + ), 
(r\E~) with their analytical dependence on the energy, and see why applying 
the Sturm-Liouville theory to these eigenfunctions yields three different RSPs. 

(iv) When two eigenfunctions have just a different normalization, they gen- 
erate the same RSP. A criterion to check whether or not two eigensolutions lead 
to different RSPs (i.e., whether or not two eigensolutions differ from a normal- 
ization factor) is provided. We shall apply the criterion to (r\E), (r\E + ), (r\E~) 
and see (as expected) that they are not a normalization of each other. 

(v) We shall compare the boundary conditions satisfied by the Gamow vec- 
tors with those satisfied by (r\E), (r\E + ), and (r\E~). More precisely, we shall 
compare the purely outgoing boundary condition with the standing-wave, "in," 
and "out" boundary conditions, and see why the purely outgoing boundary con- 
dition determines the physical content of the Gamow vectors. 

(vi) We shall discuss the time asymmetry of QED. We shall refer to this time 
asymmetry as the QED arrow of time. We shall see that the QED arrow of time is 
built into Feynman's propagator. A parallel of the QED arrow of time with that 
of Classical Electrodynamics will be made. We shall conclude that, in general, 
there exists an arrow of time at the microscopic level, and that this arrow of time 
arises from the imposition of a time-asymmetric, causal boundary condition upon 
the (time-symmetric) Schrodinger equation. Because solving the Schrodinger 
equation subject to a time-asymmetric boundary condition necessarily involves 
the construction of a propagator, we shall conclude that the quantum-mechanical 
arrow of time is built into the propagators. 

2 Boundary Conditions upon the Time-Independent 
Schrodinger Equation 

We proceed now to see how boundary conditions affect the behavior of the 
solutions of the time-independent Schrodinger equation. Rather than working 
in a general fashion, we shall use the spherical shell potential as an illustrative 
example. Generalizations to more complicated potentials are straightforward. 
Consider the spherical shell potential of height Vq, 



If we restrict ourselves to the zero angular momentum case, then the spherical 
shell Hamiltonian acts, in the radial position representation, as the following 





4 Rafael de la Madrid 

formal differential operator: 

h 2 d 2 

The time-independent Schrodinger equation (for zero angular momentum) reads 

h 2 d 2 



2mh2 +V(r))a(r;E)=Ea(r;E). (4) 

Our objective in this section is to solve Q subject to various boundary con- 
ditions, and to analyze the physical content of both the solutions and the bound- 
ary conditions. We shall study three cases: the standing-wave, the Lippmann- 
Schwinger, and the Gamow eigensolutions. 



2.1 Standing- Wave Eigenfunctions 

We first study the standing- wave eigenfunctions. To obtain them, we solve Q 
under the following boundary conditions: 

a(0;E)=0, (5a) 
a(r; E) is continuous at r = a and at r — b , (5b) 

— a(r; E) is continuous at r — a and at r = b . (5c) 
dr 

The eigensolution of (@J that satisfies l(5li|l - (f5*c)l is given by the regular solution: 

{sin(fcr) < r < a 

Ji (fc)e^ r + J 2 {k)e- i( ^ r a<r<b (6) 

J 3 (k)c ikr + Ji {k)e~ ikr b<r<oo, 

where 



The coefficients Ji(k)-J4,(k) in © are such that x( r 'iE) satisfies l|5b|l and Q5c[l. 
The expressions of Ji(k)-J^(k) can be easily calculated (they can also be found, 
for example, in pQ). 

The regular solution x( r > E) is not (5-normalized. In order to (5-normalize it, 
we need the spectral measure 

1 2m/h 2 1 . , 

Multiplying x{ r \ E) by the square root of this spectral measure yields the 6- 
normalized eigensolution of J3J that satisfies the boundary conditions (|5a(l - l(5*c| : 



(r\E) = y/g(k)x(r;k). (9) 
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The (5-normalization of (r\E) is to be understood in the following sense: 

/>oo 

/ dr(E\r)(r\E') = 6{E - E') , E, E' 6 [0, oo) , (10) 
Jo 

where (E\r) is the complex conjugate of (r\E). Note that although (r\E) is also 
defined for complex energies, we have restricted E to [0, oo), because the (Hilbert 
space) spectrum of the Hamiltonian is the positive real line. 
At infinity, the eigensolution is a linear combination of 

y/^)Mk)e- ikr , (11) 

which is an incoming spherical wave of amplitude \J g(k) ^Ti(k), and 

\^}Ja(k)e ikr , (12) 

which is an outgoing spherical wave of amplitude yj g(k) Js(k). It is easy to see 
that when E is real, l|ll|) is the complex conjugate of (|12|l . Thus, far away from 
the potential region, (r\E) is the linear combination of an incoming spherical 
wave and its complex conjugate. This behavior is very much like that of a sinu- 
soidal function - hence the name standing-wave solution for the eigenfunction 
(r\E). 

As shown in p^, the eigenfunctions 10 generate, by means of the Sturm- 
Liouvillc theory, 1 an RSP; that is, the (r\E) generate an eigenfunction expan- 
sion, a unitary operator U that diagonalizes the Hamiltonian, a direct integral 
decomposition, and an RHS 

*CWC# X . (13) 

(The explicit form of the RSP generated by the eigenfunctions © can be found 
in p.) Now, the boundary conditions (|5a|) - (|5c|) completely determine the radial 
dependence of the eigensolution of the Schrodinger equation . Essentially, the 
regular solution %(r; E) is unique up to multiplication by a function of the energy. 
Since in Quantum Mechanics the boundary conditions l|5a|) - (|5c[) are customarily 
imposed upon the Schrodinger equation, one may be tempted to conclude that 
we have found all the possible RSPs of the spherical shell potential. As we shall 
see, this is not the case: when we multiply x( r i E) by the Jost functions, we 
obtain the "in" and "out" eigensolutions, which generate two new RSPs for 
the spherical shell potential. The "in" and "out" eigensolutions (and therefore 
their associated RSPs) are determined by the boundary conditions built into the 
Lippmann-Schwingcr equation. 



2.2 Lippmann-Schwinger Eigenfunctions 

The Lippmann-Schwinger equation is usually written as 

1 Our basic reference on the Sturm-Liouville theory is |14) . Illustrative applications of 
the Sturm-Liouville theory can be found, for example, in 1 2 31151161 J7| . 
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where Hq is the free Hamiltonian. In the radial position representation, 11411 
reads 

(r|^) = (r\E) + (r| - \ V\E±) . (15) 
H> — riQ ± l£ 

The solutions (■r|£ ,± ) of this equation will be called the "in" (+) and "out" (— ) 
Lippmann-Schwinger eigenfunctions. As is well known, the integral equation 115JI 
is equivalent to the Schrodinger equation 

2^d^ + ^ (r) j (r|S >=W ±> (16) 

subject to the following boundary conditions: 

<0| JS? ± > = , (17a) 
(r|i? ± } is continuous at r = a and at r = b , (17b) 

-^-(r|i? ± ) is continuous at r — a and at r = b, (17c) 
dr 

(r|£ + ) - e"^ - S(k)e ikr as r oo , (17d) 

(r|£T) -e ifcr -5*(fc)e- ifer asr^oo, (17e) 

where S(k) is the S matrix, which is given by the quotient of the Jost functions: 

s(£,ss(t) =s!- (i8) 

The Jost functions can be written in terms of the coefficients of |JBJ as 

J+(k) = -2iJi(k); J-(k) = 2iJ 3 (k) . (19) 

The 5-normalized [in the sense of (|10|) ] Lippmann-Schwinger eigenfunctions 
can be easily obtained from H16(l - (|19|l : 

(r| £ ;±) = ^±(fc)^| i (20) 
where Q (k) are spectral measures, 

g +{k) = g-{k)= 1 - 2 -^. (21) 

7T k 

As in the standing- wave case, we are restricting E to [0, oo). From l|18l) . H2U|) and 
<|21|) it follows that the Lippmann-Schwinger eigenfunctions arc proportional to 
each other: 

(r\E+)=S(E)(r\E-). (22) 

Applying the Sturm-Liouville theory to (r|_E ± ) yields two other RSPs pQ, 
i.e., two eigenfunction expansions, two unitary operators U± that diagonalize 
the Hamiltonian, two direct integral decompositions, and two RHSs 2 

$±C«C^. (23) 



2 The RHSs <12'M are sketched in Their complete characterization is a matter of 
current investigation. 
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Therefore, even though (r\E), (r\E + ), (r\E~) all fulfill one and the same Schrodinger 
equation and have one and the same radial dependence, they generate three dif- 
ferent RSPs. Moreover, (r\E + ) and (r\E~) differ from each other by just a phase 
factor [see because ^(-B)) = 1 when E £ [0, oo). How is then possible that 

they lead to different RSPs? The answer to this question is the following: (r\E), 
(r\E + ), (r\E~) lead to different RSPs because they satisfy different, physically 
non-equivalent boundary conditions. 

In order to better understand why (r\E), (r\E + ), (r\E~) yield three different 
RSPs, we compare the Lippmann-Schwinger boundary conditions (I17a|l - I|17cll 
to the standing- wave boundary conditions (|5a|l - lfKc| . We can see first that the 
boundary conditions 117a|) - (|17c|) are the same as l|5a |) -l|5c | ). We can also see 
that in the Lippmann-Schwinger case, we have imposed an additional boundary 
condition that selects the asymptotic behavior of the eigenfunctions at infinity. 
For the "in" Lippmann-Schwinger eigenfunction, we have chosen l|17d(l . which 
means that far away from the potential region (r\E + ) is a linear combination of 
an incoming spherical wave and an outgoing spherical wave multiplied by the S 
matrix (which is a phase factor). For the "out" Lippmann-Schwinger eigenfunc- 
tion, we have chosen l|17ejl . which means that far away from the potential region 
(r\E~) is a linear combination of an outgoing spherical wave and an incoming 
spherical wave multiplied by the complex conjugate of the S matrix (which is 
also a phase factor). In the standing- wave case, we did not explicitly impose any 
boundary condition at infinity, which is tantamount to imposing the standing- 
wave asymptotic behavior. Thus (r\E), (r\E + ), (r\E~) differ from each other 
just by their asymptotic behavior at infinity. These different asymptotic behav- 
iors lead, by means of the Jost functions, to different analytical properties of 
the eigensolutions as functions of the energy when E is allowed to be complex. 
Because the Sturm-Liouville theory |14| always deals with complex energies, 
eigenfunctions with different analytical properties yield different RSPs. There- 
fore, what ultimately makes (r\E), (r\E + ), (r\E~) different is their different 
analytical behavior when E is allowed to be complex. 

An important conclusion can be drawn from the previous paragraph: bound- 
ary conditions that in the position representation select the asymptotic behavior 
read, in the energy representation, as boundary conditions that select the ana- 
lytical behavior, and vice versa. This is particularly apparent in the Lippmann- 
Schwinger equation, where the asymptotic boundary conditions (|17d|) and (|17e|) 
are built into the ±k of (|15|l . and vice versa. We then say that the analytical 
boundary conditions of the Lippmann-Schwinger equation select what is "in" 
(+ie) and what is "out" {—is) or, equivalently, that the asymptotic behaviors 
select what is "in" (I17d|) and what is "out" i|17c|l . 

The eigenfunctions (r\E), (r\E + ), (r\E~) all are proportional to the regular 
solution x{ r ', E): 



(r\E) 



(24) 



(r\E 




(25) 
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{AE ~ )= jW x(r ' ,ky (26) 

As noted above, the analytical properties of the functions that multiply x{ r \ E) 
in (|24|) - (|26[1 is what ultimately leads to different RSPs. However, it is not al- 
ways true that multiplying x( r lE) by a function of E yields an eigensolution 
that generates a different RSP. For instance, the radial solution x(r; E) and the 
eigenfunction (r\E) both lead to the same RSP. This is why we say that (r\E) is 
the ^-normalization of x(r; E). In this case, \J g(k) is just a normalization factor. 

We may then ask: given the regular solution x(r; E) and a function of the 
energy /(£"), how can we know whether f(E)x(r;E) is just a normalization 
of x( r 'i E) or leads to a different RSP? A general answer to this question is 
not known. Of course, to know the answer one can always apply the Sturm- 
Liouville theory and see if f(E)x(r;E) yields the same RSP as x{ r \E). This 
may be impractical, though. A faster method, that works at least for simple 
potentials |H15j , is to check whether 

[f(E*)]* = f(E), EeC, (27) 



or 

[f(E*)]*^f(E), EeC. (28) 

If f(E) fulfills then f(E) X (r; E) is just a normalization of x(r; E). If f(E) 
fulfills |gSJ|, then f(E)x(r;E) and x(r] E) lead to different RSPs and therefore 
have different physical content. 

To check that the criterion (|27|I - H28|) does indeed work for the spherical shell 
potential, we apply it to (r\E), (r\E + ), (r\E~). We need first to choose the 
following branch for the square root function: 

V 7 : {E e C | - 7T < arg(E) < tt} i — > {E e C | - tt/2 < arg(£) < tt/2} . (29) 

For (r\E), one can easily check that 

[q(E*)]* = q(E) , EeC. (30) 

From (J2U and (J20J, and from the criterion (27 } -{SB )l it follows that (r\E) is just 
a normalization of x(t;E). For (r|£' + ), we have that 

[g+(E*)]* = Q + (E), EEC, (31) 

but 

[J+{E*)\* = J_{E) ± J+(E) , EeC. (32) 

From (EU and (J22J), and from the criterion fZQ-fZjft it follows that (r\E+) 
is not a normalization of x(r; E) but rather has a different physical content. 
Similarly, it can be seen that the physical content of (r\E~) is not the same as 
that of (r\E). 
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The Gamow vectors are the state vectors of resonances |5lfil7l8l9ll()ll| . Like 
the standing-wave and the Lippmann-Schwinger eigensolutions, they solve the 
Schrodinger equation. At infinity, however, the Gamow eigenfunctions satisfy 
a boundary condition that is different from those satisfied by (r\E), (r\E + ), 
(r\E~): the purely outgoing boundary condition. This purely outgoing behavior 
determines the physical content of the Gamow vectors. 

There are two kinds of Gamow vectors. The first kind is the so-called decaying 
Gamow ket \z^), which is associated to a complex energy z R = Er — lTr/2 that 
lies in the lower half plane of the second sheet of the Riemann surface. The 
corresponding wave number lies in the fourth quadrant of the complex wave- 
number plane. The second kind of Gamow vector is the so-called growing Gamow 
ket |zj^ + ), which is associated to an energy zj^ = E^+U\/2 that lies in the upper 
half plane of the second sheet of the Riemann surface. The corresponding wave 
number lies in the third quadrant of the complex wave-number plane. As far as 
the time-independent Schrodinger equation is concerned, any complex number 
can be an eigenvalue of the Hamiltonian ^2]. The role of the purely outgoing 
boundary condition is to select, among all the complex energies, those that are 
to correspond to resonance energies [H]- 

In order to obtain the Gamow vectors, we solve the Schrodinger differential 
equation 

+ <r|*0 = Mr\*R) , (33) 

subject to purely outgoing boundary conditions: 
(0|«a) - 

(r|zR.) is continuous at r = a and at r = b 
— (Hzr) is continuous at r — a and at r = b 
(r|z R ) ~ e lfcRr as r — * oo , 

where 

/ 2 777 / 2 777 

kR = y~^2 ZR ' Qr= y-ffizR-Va) . (35) 

In and (|34ajHI|34d|) . (r|z R ) can denote either (r|%) or (r\z^ + ). This will 
cause no confusion, because whenever the complex energy lies in the lower half 
plane, (r|zR,) will denote (r|z^), and whenever it lies in the upper half plane, 
(r|z R ) will denote (r|zg+). 

For the spherical shell potential, l|33|l subject to the boundary conditions 
H34a|) - (|34d(> has solutions only for a denumerable set of complex energies. These 
energies come in complex conjugate pairs z n , z*, where z n = E n — iT n /2 is the 
decaying pole, and z* = E n + iP n /2 is the growing pole. The corresponding 
decaying and growing wave numbers are given by 



(34a) 
(34b) 

(34c) 

(34d) 
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In terms of the wave number k n , the nth decaying Gamow eigensohition reads 

jjj. n ) sin(fc„r) < r < a 

( 'i - = A << i flfel eiQ " r + f{fel e " iQ " r a < r < b ( 3? ) 
e ife " r b < r < oo , 

where N n is a normalization factor, 

N^ires[S(k)] k=kn , (38) 

and where 

Qn = ]/-^(Zn-V ) . (39) 

The nth growing Gamow eigensolution reads 



J* 



sin(-fc*r) 0<r<a 



M„ { f^§e^> + ^jge-«;- a < r < b (40) 
e~ lfc ™ r 6 < r < oo . 

where M„ is a normalization factor, 

A^=ires[5(fc)] fc= _ fc; = , (41) 

and where 

Qn = \/^te-Vo). (42) 
We compare now the boundary conditions (|34a|l - lj34d|l satisfied by the Gamow 



eigenfunctions with those satisfied by the standing-wave eigenfunctions [see (|5ap 
JScJ] and by the Lippmann-Schwinger eigenfunctions [sec 117al) (|17e|) ]. Clearly, 
the boundary condition that is specific to the Gamow vectors is (|34d|l . This 
boundary condition singles out the complex resonance spectrum of the Schrodinger 
equation by specifying the asymptotic behavior of the Gamow eigenfunctions. 
Moreover, the resonance spectrum selected by (|34dll coincides with the poles of 
the S matrix l(T%|l . 

To finish this section, we note that in the S'-matrix formalism, a resonance 
energy is reached by analytic continuation of S(E) from its values on the phys- 
ical spectrum to the resonance pole. In a similar vein, Gamow vectors can be 
viewed as the solutions of the analytic continuation of the Schrodinger equa- 
tion [subject to the boundary conditions (|34a|l - lj34dp ] from the energies in the 
physical spectrum to the complex resonance eigenvalue. 



3 The Arrow of Time of Quantum Electrodynamics 



We have seen how boundary conditions determine the physical content of the 
eigensolutions of the time- independent Schrodinger equation. We now turn to see 
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how boundary conditions affect the physical content of the solutions of the time- 
dependent Schrodinger equation. We shall see that time-asymmetric boundary 
conditions imposed upon the (time-symmetric) Schrodinger equation yield so- 
lutions that have an arrow of time built into them. Quantum Electrodynamics 
(QED) will be used as an illustrative example. 

In order to make a parallel with the quantum case, we recall first the essential 
features of the arrow of time of Classical Electrodynamics (CED). The Maxwell 
equations, which describe the classical electromagnetic fields, are time symmet- 
ric. The solutions of the Maxwell equations can be written as a combination 
of a retarded and an advanced solution. Experimentally, we always observe that 
light has a retarded behavior - light cannot be detected at a distance R from the 
source at any time less than R/c. In order to account for this retarded behavior, 
we select the retarded solution of the Maxwell equation and forbid the advanced 
solution. This amounts to imposing a time-asymmetric, retarded, causal bound- 
ary condition upon the (time-symmetric) Maxwell equations. Essentially, this is 
the radiation arrow of time. We stress that this arrow of time stems ultimately 
from the imposition of a causal boundary condition. 

QED is the quantum counterpart of CED, and also has an arrow of time built 
into it. The arrow of time of QED is built into the Feynman propagator. Although 
the Schrodinger equation is time symmetric, one always imposes Feynman's +ie 
prescription to construct the Feynman propagator - particles travel forward 
in time, whereas antiparticles "travel backward" in time ^3]. Thus Feynman's 
+ie prescription imposes a retarded (advanced) condition on the time evolution 
of particles (antiparticles). Essentially, this is the arrow of time of QED. In 
particular, this arrow of time sheds light onto the physical meaning of Feynman's 
+ie prescription: this prescription is just a causal boundary condition imposed 
upon the time evolution of particles and antiparticles. 

The example of QED can be generalized to any closed quantum system. 
Although the Schrodinger equation (which describes the evolution of a closed 
quantum system) is time symmetric, physical processes seem to comply with our 
basic notions of causality: cause is prior to effect. In order to make the solutions 
of the (time-symmetric) Schrodinger equation comply with causality, we impose 
causal, time-asymmetric boundary conditions upon the Schrodinger equation. 
These boundary conditions single out the causal solutions that account for the 
observed time-asymmetric phenomena. Building that (e.g., retarded) boundary 
condition into the time evolution of the quantum system involves the construc- 
tion of an (e.g., retarded) propagator. Actually, although it is well known that 
propagators have an arrow of time built into them, it is not so well empha- 
sized |11| that this implies the existence of a fundamental time asymmetry at 
the microscopic level. 

Many authors have realized the central role played by boundary conditions 
in the description of time asymmetry. For example, Ritz |19j thought that the 
time asymmetry (often called irreversibility) of statistical mechanics arises from 
boundary conditions (in contrast to Einstein |19| . who thought that irreversibil- 
ity comes from averaging over a large number of systems, that is, for Einstein ir- 
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reversibility emerges from probability and statistics). For Peierls [201 irreversibil- 
ity also arises as a consequence of boundary conditions. To show this, Peierls 
rephrases Boltzmann's Stosszahl-Ansatz, which is the origin of the irreversibility 
of the so-called "Lorentz g 3jS, cLS cL boundary condition [201 - Some authors such 
as Penrose [21] or Gell-Mann and Hartle [22] have used boundary conditions 
as a possible explanation of time asymmetry, and even of time symmetry |22) . 
Other authors such as Preskill [22| also use boundary conditions to explain the 
time asymmetry (irreversibility) of quantum statistical mechanics: essentially, 
the non-decreasing behavior of the entropy can be understood as stemming 
from the assumption that system and environment are initially uncorrelated, 
i.e., from the assumption that the initial system-environment state is separable 
(unentangled) [23j- 

4 Conclusions 

We have seen how important boundary conditions are in obtaining the solutions 
of the time-independent Schrodinger equation that fit a given physical situation. 
Essentially, the asymptotic behavior of the eigensolution determines its physical 
content. We have analyzed and compared the standing-wave, "in," "out," and 
purely outgoing boundary conditions. We have seen that the standing-wave, 
"in," and "out" boundary conditions yield three physically different RSPs. The 
purely outgoing boundary condition selects the resonance spectrum. 

We have discussed the time-asymmetry of QED. Essentially, Feynman's +ie 
prescription, which is used to construct the Feynman propagator, encodes the 
time asymmetry of QED. We have also concluded that in general, the time 
evolution of a closed quantum system has a time asymmetry built into the prop- 
agators. 
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